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LARGE  EDDY  SIMULATION  OF 
THE  FLOW  OVER 

A  LOW  REYNOLDS-NUMBER  AIRFOIL 
Final  report 

Ugo  Piomelli 

Dept,  of  Mechanical  Engineering 
University  of  Maryland 
College  Park 

During  the  grant  period,  April  1,  1997-December  31,  1998,  the  following  tasks  were 
accomplished: 

1.  Model  development 

An  a  priori  study  of  subgrid-scale  models  for  the  unclosed  terms  in  the  energy  equation 
has  been  carried  out  using  the  flow  field  obtained  from  the  direct  simulation  of  homogeneous 
isotropic  turbulence.  Scale-similar  models  involve  multiple  filtering  operations  to  identify  the 
smallest  resolved  scales  that  have  been  shown  to  be  the  most  active  in  the  interaction  with  the 
unresolved  subgrid  scales  (SGS).  These  models  were  found  to  give  more  accurate  prediction 
of  the  SGS  stresses  and  heat  fluxes  than  eddy- viscosity  and  eddy-diffusivity  models,  as  well 
as  improve  prediction  of  the  SGS  turbulent  diffusion,  SGS  viscous  dissipation,  and  SGS 
viscous  diffusion.  The  turbulent  diffusion  term  has,  in  the  past,  been  neglected,  or  modeled 
heuristically.  The  present  work  examined  this  term  rationally  for  the  first  time.  The  model 
proposed  is  inexpensive  (since  it  does  not  require  the  calculation  of  additional  terms,  but  is 
based  on  the  model  fo  the  SGS  stresses  themselves),  and  gives  good  results  in  the  a  priori 
tests  carried  out  so  far. 

This  work  has  resulted  in  a  paper  accepted  at  the  3rd  ASME/JSME  Joint  Fluids  En¬ 
gineering  Conference,  which  was  subsequently  expanded  and  submitted  to  Theoretical  and 
Computational  Fluid  Dynamics.  The  two  papers  are  enclosed. 

2.  Boundary-layer  calculations 

Time-developing  turbulent  boundary  layers  over  an  isothermal  flat  plate  at  free-stream  Mach 
numbers  of  0.3  and  0.7  have  been  computed  using  an  explicit  finite-difference  method  on 
structured  multi-block  grids.  The  size  of  each  block  is  adjusted  depending  on  the  dimension 
of  the  largest  structures  present  locally  in  the  flow.  This  results  in  substantial  savings 
of  memory  and  CPU  time,  compared  to  standard  single-block  methods.  In  the  calculations 
presented  the  near-wall  region  is  computed  using  a  domain  with  a  spanwise  length  Lj"  =  820, 
which  is  sufficient  to  contain  several  streaks.  This  grid  block  is  repeated  periodically  in  the 
spanwise  direction.  The  outer  layer,  which  contains  larger  structures,  is  computed  using  a 
domain  that  is  twice  as  wide  (L+  =  1640).  Although  the  flow  at  the  interface  between  the 
blocks  has  a  periodicity  length  determined  by  the  inner-layer  block,  within  a  few  grid  points 
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longer  wavelengths  are  generated  by  the  non-linear  interactions.  The  velocity  statistics  and 
rms  intensities  compare  well  with  single-block  calculations  that  use  substantially  more  grid 
points. 

The  advantage  of  this  method  is  that,  for  flows  that  are  homogeneous  in  one  direction,  it 
makes  the  cost  of  the  calculation  nearly  independent  of  the  Reynolds  number,  without  requir¬ 
ing  additional  modeling  of  the  wall  layer.  Flows  of  this  type  include  airfoils,  axisymmetric 
combustors,  and  other  two-dimensional  or  axisymmetric  conflgurations. 

This  work  has  resulted  in  a  paper  that  has  been  submitted  to  the  Journal  of  Computa¬ 
tional  Physics,  and  that  will  also  be  presented  at  the  Second  Air  Force  International  Direct 
and  Large-eddy  Simulation  Conference.  The  paper  is  enclosed. 


Proceedings  of  the 
3rd  ASME/JSME  Joint  Fluids  Engineering  Conference 
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ABSTRACT 

An  a  priori  study  of  subgrid>scale  models  for  the  unclosed 
terms  in  the  energy  equation  is  carried  out  using  the  flow  field 
obtained  from  the  direct  simulation  of  homogeneous  isotropic 
turbulence.  Scale-similar  models  involve  multiple  filtering  op¬ 
erations  to  identify  the  smallest  resolved  scales  and  have  been 
shown  to  be  the  most  active  in  the  interaction  with  the  unre¬ 
solved  subgrid  scales  (SGS).  In  the  present  study  these  models 
are  found  to  give  more  accurate  prediction  of  the  SGS  stresses 
and  heat  fluxes  than  eddy- viscosity  and  eddy-diffusivity  mod¬ 
els.  as  well  as  improve  prediction  of  the  SGS  turbulent  diffusion, 
SGS  viscous  dissipation,  and  SGS  viscous  diffusion. 

1  Introduction 

Large-eddy  simulation  (LES)  is  a  technique  intermediate  be¬ 
tween  the  direct  simulation  (DNS)  of  turbulent  flows  and  the  so¬ 
lution  of  the  Reynolds-averaged  equations.  In  LES  the  contribu¬ 
tion  of  the  large,  energy-carrying  structures  to  momentum  and 
energy  transfer  is  computed  accurately,  and  only  the  effect  of  the 
smallest  scales  of  turbulence  is  modeled.  Since  the  small  scales 
tend  to  be  more  homogeneous  and  universal,  and  less  affected  by 
the  boundary  conditions  than  the  large  ones,  there  is  hope  that 
their  models  can  be  simpler  and  require  fewer  adjustments  when 
applied  to  different  flows  than  similar  models  for  the  Reynolds- 
averaged  Navier-Stokes  equations. 

While  a  substantial  amount  of  research  has  been  carried  out 
into  the  modeling  aspects  and  requirements  for  incompressible 
flows,  the  applications  of  large-eddy  simulation  to  compress- 
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ible  flows  have  been  significantly  fewer.  One  of  the  reasons 
for  the  comparatively  small  number  of  calculations  of  compress¬ 
ible  flows  is  undoubtedly  the  additional  complexity  introduced 
by  the  need  to  solve  an  energy  equation,  which  introduces  ex¬ 
tra  unclosed  terms.  In  addition  to  the  subgrid  scale  stresses  that 
must  be  modeled  in  incompressible  flows  as  well,  several  other 
unclosed  terms  appear  in  the  filtered  equations  for  compressible 
flows.  Furthermore,  the  form  of  the  unclosed  terms  depends  on 
the  form  of  the  energy  equation  chosen  (internal  or  total  energy, 
total  energy  of  the  resolved  field  or  enthalpy). 

Early  applications  of  LES  to  compressible  flows  used  a 
transport  equation  for  the  internal  energy  per  unit  mass,  e  (Moin 
et  ai  1991,  El-Hady  et  al  1994)  or  for  the  enthalpy  per  unit 
mass,  h  (Speziale  et  al.  1988,  Erlebacher  et  al.  1992).  In  these 
equations,  the  SGS  heat  flux  was  modeled  in  a  manner  similar  to 
that  used  for  the  SGS  stresses,  and  the  remaining  terms  (the  SGS 
pressure-dilatation  flj//,  and  the  SGS  contribution  to  the  viscous 
dissipation,  Ev)  were  neglected. 

Vreman  et  al.  (1995b)  performed  a  priori  tests  using  DNS 
data  obtained  from  the  calculation  of  a  mixing  layer  at  Mach 
numbers  in  the  range  0.2-0.6  to  establish  the  validity  of  these 
assumptions.  They  found  that  the  SGS  pressure-dilatation  and 
SGS  viscous  dissipation  are  of  the  same  order  as  the  divergence 
of  the  SGS  heat  flux  Qj\  and  that  modeling  ty  improves  the  re¬ 
sults,  especially  at  moderate  or  high  Mach  numbers. 

Vreman  etai  {\ 995a,  1 995b)  proposed  the  use  of  a  transport 
equation  for  the  total  energy  of  the  filtered  field,  rather  than  either 
the  enthalpy  or  the  internal  energy  equations;  the  same  unclosed 
terms  that  appear  in  this  equation  are  also  present  in  the  internal 
energy  and  enthalpy  equations.  This  equation  was  also  used  by 
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Normand  and  Lesieur  (1992),  who  neglected  both  Tldu  and  Ev,.. 

Very  few  calculations  have  been  carried  out  using  the  trans¬ 
port  equation  for  the  total  energy,  despite  the  desirable  feature 
that  it  is  a  conserved  quantity,  and  that  all  the  SGS  terms  in  this 
equation  can  be  cast  in  conservative  form.  This  equation  has 
a  completely  different  set  of  unclosed  terms,  whose  modeling 
is  not  very  advanced  yet.  Knight  et  ai  (1998)  performed  the 
LES  of  isotropic  homogeneous  turbulence  on  unstructured  grids 
and  compared  the  results  obtained  with  the  Smagorinsky  (1963) 
model  with  those  obtained  when  the  energy  dissipation  was  pro¬ 
vided  only  by  the  dissipation  inherent  in  the  numerical  algorithm. 
They  modeled  the  SGS  heat  flux  and  an  SGS  turbulent  diffusion 
term,  and  neglected  the  SGS  viscous  diffusion. 

In  this  paper,  the  flow  field  from  a  DNS  of  homogeneous 
isotropic  turbulence  is  used  to  compute  the  terms  in  the  energy 
equations,  and  evaluate  possible  models  for  their  parameteriza¬ 
tion.  The  work  will  be  focused  mainly  in  the  total  energy  equa¬ 
tion.  both  because  of  the  lack  of  previous  studies  of  the  terms  that 
appear  in  it,  and  because  of  the  desirability  of  solving  a  transport 
equation  for  a  conserved  quantity. 

In  the  following  section,  the  governing  equations  are  pre¬ 
sented,  the  terms  that  require  closure  are  defined,  and  the  DNS 
database  used  for  the  a  priori  tests  is  described.  In  Sections  3 
and  4  several  models  for  the  unclosed  terms  are  presented  and 
tested.  Finally,  some  conclusions  are  drawn  in  Section  5. 


tions  of  conservation  of  mass,  momentum  and  energy.  In  com¬ 
pressible  flows,  in  addition  to  the  mass  and  momentum  equa¬ 
tions,  one  can  choose  between  solving  an  equation  for  the  inter¬ 
nal  energy,  enthalpy  or  total  energy.  The  filtered  equations  of 
motion,  then,  can  be  put  in  the  form: 
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2  Problem  formulation 
2.1  Governing  equations 

To  separate  the  large  from  the  small  scales,  LES  is  based 
on  the  definition  of  a  filtering  operation:  a  filtered  (or  resolved, 
or  large-scale)  variable,  denoted  by  an  overbar,  is  defined  as 
(Leonard  1974) 

7(x)  =  [  f{x')G{\,x'\A)d\',  (1) 

Jo 

where  D  is  the  entire  domain  and  G  is  the  filter  function,  and  A  is 
the  filter  width,  z.e.,  the  wavelength  of  the  smallest  scale  retained 
by  the  filtering  operation.  The  filter  function  determines  the  size 
and  structure  of  the  small  scales. 

In  compressible  flows,  it  is  convenient  to  use  Favre-filtering 
(Favre  1965a,  1965b)  to  avoid  the  introduction  of  subgrid-scale 
terms  in  the  equation  of  conservation  of  mass.  A  Favre-filtered 
variable  is  defined  as: 

/  =  ^/P-  (2) 

To  obtain  the  equations  governing  the  motion  of  the  resolved  ed¬ 
dies,  the  Favre-filtering  operation  must  be  applied  to  the  equa¬ 


Here,  uj  is  the  velocity  in  the  j  direction,  p  is  the  density,  p  the 
pressure  and  T  the  temperature,  t  —  CyT  is  the  internal  energy 
per  unit  mass,  /i  =  e-f  p/p  is  the  enthalpy  per  unit  mass.  E  = 
e-l-  UiUi/2  is  the  total  energy  per  unit  mass,  and  the  diffusive 
fluxes  are  given  by 


where  p  is  the  molecular  viscosity,  and  k  is  the  thermal  conduc¬ 
tivity  corresponding  to  the  filtered  temperature  7.  The  effect  of 
the  subgrid  scales  appears  through  the  SGS  stresses  X/y,  the  SGS 
heat  flux  Qj,  the  SGS  pressure-dilatation  flj//,  the  SGS  contri¬ 
bution  to  the  viscous  dissipation,  e^,  the  SGS  turbulent  diffusion 
djj /dxj,  and  the  SGS  contribution  to  viscous  diffusion,  d(Dj/dxj; 
these  quantities  are  defined  as: 


=  p{uiuj-UiUj) 

(9) 

=  p{ujT-Ujf^ 

(10) 
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(11) 
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Figure  1 .  Energy  spectra  at  //T|  =  1 , 3  and  4.  The  squares  correspond 
to  the  filter-widths  used  in  the  a  priori  tests. 

The  equation  of  state  has  been  used  to  express  pressure-gradient 
and  pressure-diffusion  correlations  in  terms  of  Qj  and  Uju.  It  is 
also  assumed  here  that 

(15) 

and  an  equivalent  equality  involving  the  thermal  conductivity  ap¬ 
plies.  Vreman  et  ai  (1995b)  performed  a  priori  tests  using  DNS 
data  obtained  from  the  calculation  of  a  mixing  layer  at  Mach 
numbers  in  the  range  0.2-0. 6,  and  concluded  that  neglecting  the 
nonlinearities  of  the  diffusion  terms  in  the  momentum  and  en¬ 
ergy  equations  is  acceptable. 

2.2  A  priori  tests 

One  method  to  evaluate  the  performance  of  models  for  LES 
or  RANS  calculations  is  the  a  priori  test,  in  which  the  velocity 
fields  obtained  from  a  direct  simulation  are  filtered  to  yield  the 
exact  SGS  terms,  and  the  filtered  quantities  are  used  in  a  model¬ 
ing  ansatz  to  evaluate  the  accuracy  of  the  parameterization.  The 
database  used  in  this  study  was  obtained  from  the  calculation  of 
homogeneous  isotropic  turbulence  decay. 

The  Navier-Stokes  equations  were  integrated  in  time  using  a 
fourth-order  order  Runge-Kutta  method.  The  spatial  derivatives 
were  computed  using  an  eighth-order  accurate  central  finite- 
difference  scheme.  The  simulations  were  performed  on  grids 
with  256'^  points,  so  that  a  large  range  of  scales  is  found  in  the 
energy  spectrum.  The  computational  domain  is  a  periodic  box 
with  length  27t  in  each  dimension.  The  fluctuating  fields  were 
initialized  as  in  Martin  and  Candler  (1996)  and  the  DNS  re¬ 
sults  were  validated  by  comparison  with  the  Martin  and  Candler 
(1996)  simulations. 

The  calculation  was  performed  at  a  Reynolds  number  Re\  — 
u^X/v  =  35,  where  X  is  the  Taylor  micro-scale  and  u'  is  the  turbu¬ 


lence  intensity,  and  at  a  turbulent  Mach  number  M;  =  qja  =  0.52, 
where  =  UiUi  and  a  is  the  speed  of  sound.  Since  the  dilata- 
tional  field  is  initially  zero,  the  flow  is  allowed  to  evolve  for  one 
dimensionless  time  unit,  T/ =  X/w'. 

The  subgrid  scale  quantities  were  then  evaluated.  The  DNS 
fields  were  filtered  using  a  top-hat  filter 

+  (16) 

with  varying  filter  widths  A  =  nA,  where  A  is  the  grid  size  and  n  = 
2,  4,  6,  8  and  10.  Figure  1  shows  the  energy  spectrum  including 
the  location  of  the  filter  cutoffs.  All  the  filter- widths  tested  lie 
in  the  decaying  region  of  the  spectrum.  Most  of  the  results  will 
be  shown  for  a  filter-width  A  =  8A,  at  the  edge  of  the  inertial 
range  of  the  spectrum.  With  this  filter  width  approximately  1 17r 
of  the  energy  resides  in  the  subgrid-scales,  a  value  representative 
of  actual  LES.  Two  quantities  are  used  to  evaluate  the  accuracy 
of  a  model:  the  correlation  coefficient  of  the  modeled  term  with 
the  exact  one,  defined  as 

(/model/pNs)  (17) 

and  the  L2— norm  of  the  modeled  and  exact  terms. 


3  Models  for  the  momentum  equation 

The  modeling  of  the  SGS  stresses  has  received  compar¬ 
atively  more  attention  than  any  of  the  other  unclosed  terms 
in  compressible  flows.  Yoshizawa  (1986)  proposed  an  eddy- 
viscosity  model  for  weakly  compressible  turbulent  flows  us¬ 
ing  a  multi-scale  direct-interaction  approximation  method.  The 
anisotropic  part  of  the  SGS  stresses  is  parameterized  using  the 
Smagorinsky  (1963)  model,  while  the  SGS  energy  is 
eled  separately: 

=  -C22A‘p|5|  =  C;a„  (18) 

Ttt  =  C,2pA^|5p  =  C/a  (19) 

with  Cv  =  0.16  and  C/ =  0.09. 

Speziale  et  ai  (1988)  proposed  the  addition  of  a  scale- 
similar  part  to  the  eddy-viscosity  model  of  Yoshizawa.  Scale- 
similar  models  are  based  on  the  assumption  that  the  most  active 
subgrid  scales  are  those  closer  to  the  cutoff  wavenumber,  and  that 
the  scales  with  which  they  interact  most  are  those  right  above 
the  cutoff  (Bardina  et  ai,  1980).  The  mixed  model  proposed  by 


3 


i 
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Speziale  et  al  (1988),  and  used  by  Erlebacher  et  al  (1992)  and 
Zang  etal  (1992)  was  given  by 

§•  5' 

T’ij  -  +  Aij  -  -^Akk  (20) 

Xick  =  Qa'^Auc,  (21) 


where  Aij  =  p{uiUj  —  UiUj).  Erlebacher  et  at.  (1992)  tested  the 
constant  coefficient  model  a  priori  and  by  comparing  DNS  and 
LES  results  of  compressible  isotropic  turbulence  and  found  good 
agreement  in  the  dilatational  statistics  of  the  flow,  as  well  as  high 
correlation  between  the  exact  and  the  modeled  stresses.  Zang  et 
ai  (1992)  compared  the  DNS  and  LES  results  of  isotropic  tur¬ 
bulence  with  various  initial  ratios  of  compressible  to  total  kinetic 
energy.  They  obtained  good  agreement  for  the  evolution  of  quan¬ 
tities  such  as  compressible  kinetic  energy  and  fluctuations  of  the 
thermodynamic  variables. 

Moin  et  a/.  (1991)  proposed  a  modification  of  the  eddy- 
viscosity  model  (18-19)  in  which  the  two  model  coefficients 
were  determined  dynamically,  rather  than  input  a  priori,  using 
the  identity  (Germano  1992)  Lij  =  Tij  -  f/},  which  relates  the 
“resolved  turbulent  stresses”, 

Lij  =  (puipuj/p^-puipuj/p,  (22) 


the  subgrid-scale  stresses  1,7  and  the  subtest  stresses  Tij  = 
p  w/Wj  -  p  UjUj,  where  /  =  p//p ,  and  the  hat  represents  the 
application  of  the  test  filter  G,  of  characteristic  width  A  =  2A. 
Moin  et  al.  (1991)  determined  the  model  coefficients  by  sub¬ 
stituting  the  models  (18-19)  into  (22)  and  contracting  with  5,7; 
in  this 'work  the  contraction  proposed  by  Lilly  (1992)  to  min¬ 
imize  the  error  in  a  least-squares  sense  will  be  used  instead. 
Accordingly,  the  two  model  coefficients  for  the  Dynamic  Eddy- 
Viscosity  model  (denoted  hereafter  by  the  acronym  DEV)  will  be 
given  by 

r  _  1  {-^im^nn)  ^  _  i^k) 

^  3  {MM  ’  '  -  (p  _  a)  V  ' 

where  (i/^  =  -2Ay  |5|(57  -  bijSkk/^).  Mij  =  p,7  -  (J  = 

2A^p  |5|^,  and  the  brackets  (•)  denote  averaging  over  the  compu¬ 
tational  volume.  Dynamic  model  adjustment  can  be  also  applied 
to  the  mixed  model  (20-21),  to  yield  the  Dynamic  Mixed  model 
(DMM) 


C  = 


0  = 


{LijM,j)  -  {NjjMjj)  _  1  (4,m,M„„)  +  (/V„,„,H.,,)  24x 

{MikM,k)  3  {M,kMik)  ■ 

(P-5>  ’  ^ 


Figure  2.  A  priori  comparison  of  the  normal  SGS  stresses  Ti  i . - - 

DEV;  - — DMM; - DMM-1;  A  DNS.  (a)  Correlation  coefficient:  (b) 

fluctuations. 


with  Bij  =  p  {uiUj  —  UiUj),  and  Nij  =  Bij  -  Ai).  One  advantage 
of  mixed  models  is  that  they  allow  one  to  model  the  trace  of  the 
SGS  stresses  without  requiring  a  separate  term  of  the  form  (19). 
A  one-coefficient  dynamic  mixed  model  (DMM-1 )  would  be  of 
the  form 

1,7  =  Caij-{-Aij,  (26) 

with 

^  {UjMij)-\NijM,j) 

(MikMik) 

Figures  2  and  3  compare  the  diagonal  and  off-diagonal  com¬ 
ponents  of  the  SGS  stress  tensor  predicted  by  the  various  mod¬ 
els.  Consistent  with  the  results  of  previous  investigators  eddy- 
viscosity  models  are  not  to  be  able  to  predict  the  rms  of  the  SGS 
stresses  very  accurately.  All  models  have  high  correlation  with 
the  DNS  data,  although,  for  the  DEV  model,  that  is  due  to  the 
trace  of  The  eddy- viscosity  prediction  of  the  off-diagonal 
terms  (Fig.  3)  has,  in  fact,  a  much  lower  correlation  coefficient. 
In  general,  the  one-coefficient  mixed  model  (26-27)  appears  to 
be  the  most  accurate  among  those  tested.  Its  correlation  with  the 
exact  SGS  stresses  is  always  greater  than' 0.8,  and  the  prediction 
of  the  rms  is  consistently  more  accurate  than  that  of  the  eddy- 
viscosity  model  (and  is  also  more  accurate  than  that  obtained 
with  the  two-coefficient  mixed  model,  in  which  the  SGS  energy 
is  modeled  separately). 

The  coefficient  Q-  remained  nearly  constant  at  a  value  of 
0.15  throughout  the  calculation,  consistent  with  the  theoretical 
arguments  (Yoshizawa  1986).  The  coefficient  of  the  SGS  en¬ 
ergy,  C/,  on  the  other  hand,  has  a  value  three  times  higher  than 
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Figure  3.  A  priori  comparison  of  the  off-diagonal  SGS  stresses  X\2> 

— —  DEV; - DMM  and  DMM-1 ;  A  DNS.  (a)  Correlation  coefficient; 

(b)  rms  fluctuations. 

predicted  by  the  theory,  consistent  with  the  results  of  Moin  et 
al  (1991). 


Figure  4.  Comparison  of  unclosed  terms  in  the  energy  equations,  (a) 
Terms  in  the  internal  energy  or  enthalpy  equations;  (b)  total  energy  equa- 

IjQP,  - Divergence  of  the  SGS  heat  flux,  dQjIdxf,  — • —  SGS 

viscous  dissipation  Ey; - pressure  dilatation  11^//; - -SGS  turbu¬ 

lent  diffusion  djj Idxy,  — A—  SGS  viscous  diffusion  32)/  jdXf. 

A  mixed  model  of  the  form 


4  Models  for  the  energy  equations 

A  comparison  of  the  magnitude  of  the  unclosed  terms  in  the 
three  forms  of  the  energy  equation  (5),  (6)  and  (7)  is  shown 
in  Fig.  4.  Unlike  in  the  mixing  layer  studied  by  Vreman  ei 
ai  (1995b),  in  this  flow  the  pressure  dilatation  Uju  is  negligi¬ 
ble,  and  the  viscous  dissipation  Ev  is  less  than  one-tenth  of  the 
divergence  of  the  SGS  heat  flux.  Thus,  the  only  term  that  re¬ 
quires  modeling  in  the  internal  energy  or  enthalpy  equations  is 
QjAn  the  total  energy  equation,  on  the  other  hand,  the  SGS  tur¬ 
bulent  diffusion  djj/dxj  is  significant.  In  the  following,  several 
models  for  the  more  significant  terms  will  be  examined. 

The  most  important  term  to  be  closed  (Fig.  4)  is  the  diver¬ 
gence  of  the  SGS  heat  flux  (10).  The  simplest  approach  is  to  use 
an  eddy-diffusivity  model  of  the  form: 


^  pvr  df  A^p\S\  df 

^  Prr  dxj  Prj  dxj' 


(28) 


where  C  is  the  eddy-viscosity  coefficient  in  (23).  The  turbulent 
Prandtl  number  Prj  can  be  fixed,  or  adjusted  dynamically  ac¬ 
cording  to  Prr  =  C{TkTk) I where 


?(,  =  (p«;  pT/p  )  -  PM;  pr/p , 

(29) 

x,  ?f  ~  af 

+Ap|51^.  , 

(30) 

,A^p|5|  dT 
Prj  dxj 


-F 


(31) 


was  proposed  by  Speziale  et  ai  (1988).  The  model  coefficient  C 
is  given  by  (24);  Prj  can  again  be  assigned  a  priori  or  adjusted 
dynamically  according  to 


PrT=C 


{TkTk) 


(32) 


with 


yj  =  9  ^jT-ujT]-p[ujT-UjT\.  (33) 


In  Figure  5  the  models  described  above  are  compared.  The 
eddy-diffusivity  models  have  poor  correlation  with  the  DNS  data, 
as  is  the  case  with  these  types  of  model.  The  constant-Pr7  case 
(the  value  used  was  0.7,  following  Zang  et  ai  1992)  gives  a  very 
low  value  of  the  rms  of  the  modeled  Qj.  A  lower  value,  Prj  — 
0.4,  as  used  by  Speziale  et  ai  (1988)  would,  however,  give  rms 
fluctuations  nearly  identical  to  those  predicted  using  the  dynamic 
procedure.  The  mixed  model  gives  the  best  correlation  with  the 
data  (again  around  0.8),  but  over-predicts  the  rms\  at  early  times 
the  mixed  model  gives  results  in  better  agreement  with  the  data 
than  the  eddy-diffusivity  one,  but  for  tjXt  >  3  the  latter  gives 
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Figure  5.  Coefficient,  correlation  and  rms  of  the  model  for  the  SGS  heat 
flux  Qj.  (a)  Turbulent  Prandtl  number,  Prj\  (b)  correlation  coefficient; 

(c)  rms.  — —  Eddy-diffusivity  model,  fixed  Prandtl  number; - eddy- 

diffusivity  model,  variable  Prandtl  number; - mixed  model;  A  DNS. 


more  accurate  results.  A  limitation  of  this  study  is  the  fact  that 
the  filter  width  is  already  in  the  decaying  region  of  the  spectrum, 
a  situation  that  has  been  shown  to  degrade  the  accuracy  of  the 
dynamic  procedure  (Lund  and  Meneveau  1997). 

The  other  term  that  can  be  significant  in  the  enthalpy  or  in¬ 
ternal  energy  equations  is  the  viscous  dissipation  Ey.  Vreman  et 
ai  (1995b)  proposed  three  models  for  this  term: 


II 

(34) 

el"’  =  Ct2pq^ /Is., 

(35) 

eP’  =  Ce3p^^A, 

q‘  '^  UkUk-UkUk- 

(36) 

The  first  is  a  scale-similar  model;  the  second  and  third  use  di¬ 
mensional  analysis  to  represent  the  SGS  dissipation  as  the  ratio 
between  the  cube  of  the  SGS  velocity  scale,  q,  and  the  length 
scale,  and  assign  the  velocity  scale  using  either  the  Yoshizawa 
(1986)  model  (19)  or  the  scale-similar  model.  For  consistency, 
each  of  the  last  two  models  should  be  coupled  with  the  corre¬ 
sponding  model  for  Xkk-  Based  on  their  DNS  data,  Vreman  et 
al.  (1995b)  fixed  the  values  of  the  coefficients  that  give  the  cor¬ 
rect  magnitude  for  this  term  and  obtained:  Cgi  =  8,  Ce2  =  1 .6  and 
Ce3  =  0.6.  Alternatively,  the  dynamic  procedure  can  be  used  to 


0 


Figure  6.  Coefficient,  correlation  and  rms  of  the  model  for  the  viscous 
dissipation  Ev  (a)  Model  coefficient;  (b)  correlation  coefficient:  (c)  rms. 

- Scale  similar  (34); - Dynamic  (35): - Dynamic  (36): 

A  DNS. 

give 

-  (37) 

where 


£<'>  =  Ce, 

5/7^/7  j  ; 

(38) 

\ 

=Ce2p^/A. 

/ 

^~A-|S|^ 

(39) 

^  ^UiUj  —  UjUj. 

(40) 

Model  coefficients  obtained  frdm  the  dynamic  procedure  in 
this  form  (in  which  there  is  no  contraction)  can  become  ill- 
conditioned,  since  the  two  terms  in  the  denominator  may  be  ap¬ 
proximately  equal,  giving  spuriously  high  values  of  the  denom¬ 
inator.  This  behavior  was  observed  in  model  (34),  in  which  ac¬ 
ceptable  results  were  obtained  only  if  Qi  was  constrained  to  be 
positive,  and  model  (35).  The  model  given  by  (36),  on  the  other 
hand,  was  well  behaved. 

Figure  6  compares  the  predictions  of  the  three  models.  The 
values  of  the  coefficients  obtained  from  the  present  a  priori  test 
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1.5 


Figure  7.  Correlation  and  mis  of  the  model  for  the  turbulent  diffusion  Jj. 

(a)  correlation  coefficient;  (b)  rms.  - —  Knight  et  al  (1998)  model; 

A  DNS. 


are  lower  than  those  obtained  in  the  mixing  layer  by  Vreman 
et  al.  (1995b).  The  pure  scale-similar  model  (34)  and  the  model 
(36),  which  also  uses  scale  similarity  to  supply  the  velocity  scale, 
give  the  best  correlation  and  nearly  correct  rms  amplitudes.  The 
rms  predicted  by  the  model  (35)  is  two  orders  of  magnitude  larger 
than  the  others,  and  cannot  be  seen  in  the  plot.  In  this  flow  the 
coefficients  obtained  from  the  mixing  layer  data  would  yield  high 
values  of  the  modeled  rms,  indicating  some  lack  of  universality 
for  the  modeling  of  this  term. 

The  two  terms  in  the  total  energy  equation  that  require  mod¬ 
eling  are  the  SGS  turbulent  diffusion  dJj/dxj  and  the  SGS  vis¬ 
cous  diffusion  dT>j/dxj.  The  only  calculation  that  attempted  to 
model  the  former  was  that  by  Knight  et  al.  (1998).  They  argued 
that  Ui  Ui  and  proposed  a  model  of  the  form 


Figure  8.  Coefficient,  correlation  and  rms  of  the  model  for  the  viscous 
diffusion  T>j.  (a)  Model  coefficient;  (b)  correlation  coefficient;  (c)  rms. 
- Scale-similar  model;  A  DNS. 


model  for  this  term  has  been  proposed  in  the  literature  to  date. 
One  possibility  is  to  parameterize  it  using  a  scale-similar  model 
of  the  form 


(Dj  —  C[){(5ijUi 


(42) 


_  in  which  the  coefficient  can  be  obtained  from 

JjC:tUkZjk.  (41) 


This  model  is  compared  in  Fig.  7  with  the  DNS  data;  Zjk  was 
obtained  from  the  DMM-1  model  (41).  The  model  has  a  high 
correlation  with  the  data  (of  order  0.8),  and  the  rms  also  matches 
the  data  well.  It  should  be  noticed,  however,  that  the  model  is 
built  upon  the  prediction  of  the  SGS  stresses  by  DMM-1,  which 
over-predicts  the  rms  of  the  normal  stresses  by  30%,  that  of  the 
off-diagonal  ones  by  about  50%.  It  appears  that  the  modeling 
assumption  by  itself  might  underestimate  the  diffusion,  an  error 
that  is  compensated  by  one  of  opposite  sign  in  the  SGS  stress 
model.  The  high  correlation,  however,  indicates  that  addition  of 
a  model  coefficient,  perhaps  adjusted  dynamically,  may  be  bene¬ 
ficial. 

The  SGS  viscous  diffusion  d(Dj/dxj  is  the  smallest  of  the 
terms  in  the  total  energy  equation.  It  is  5%  of  the  divergence  of 
Qj  SLit/Xi  =  1,  but  increases  to  about  10%  at  the  final  time.  No 


Cz)  = 


_  pg/  / P»i 

p"  p‘ 


M) 


(43) 


where 


IKl=  (oikUk-0{kU}\  -  (oikUk-OikUk 


(44) 


As  can  be  seen  from  Fig.  8,  however,  this  approach  gives  a  fairly 
poor  correlation,  and  fair  agreement  for  the  prediction  of  the  rms 
intensities.  This  error  may,  however,  be  tolerable  given  the  small 
contribution  that  this  term  gives  to  the  energy  budget. 
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5  Conclusions 

Several  mixed  and  eddy- viscosity  models  for  the  momen¬ 
tum  and  energy  equations  have  been  tested.  The  velocity,  pres¬ 
sure,  density  and  temperature  fields  obtained  from  the  DNS  of 
homogeneous  isotropic  turbulence  at  Rex  =  35,  Mt  —  0.52  were 
filtered  and  the  unclosed  terms  in  the  momentum,  internal  energy 
and  total  energy  equations  were  computed. 

In  the  momentum  equation,  mixed  models  were  found  to 
give  better  prediction,  in  terms  of  both  correlation  and  rms  ampli¬ 
tude,  than  the  pure  eddy-viscosity  models.  The  dynamic  adjust¬ 
ment  of  the  model  coefficient  was  beneficial,  as  already  observed 
by  Moin  er  a/.  (1991). 

In  the  internal  energy  and  enthalpy  equations,  in  this  flow, 
only  the  divergence  of  the  SGS  heat  flux  was  significant;  the  SGS 
pressure  dilatation  Wdu  and  viscous  dissipation  Ev,  which  were 
significant  in  the  mixing  layer  studied  by  Vreman  et  al  (1995b), 
were  found  to  be  negligible  here.  Once  again,  mixed  dynamic 
models  gave  the  most  accurate  results.  In  particular,  the  turbulent 
Prandtl  number  obtained  dynamically  was  somewhat  lower  than 
the  value  of  0.7  often  assigned  a  priori. 

In  the  total  energy  equation  two  additional  terms  are  present, 
one  of  which,  the  turbulent  diffusion  dJj/dxj  is  significant.  The 
model  proposed  by  Knight  ef  al.  (1998),  which  parameterizes 
the  turbulent  diffusion  in  terms  of  the  SGS  stresses,  correlates 
well  with  the  actual  SGS  stresses,  and  predicts  the  correct  rms 
amplitude,  A  mixed  model  for  the  SGS  turbulent  diffusion  has 
also  been  proposed  and  tested,  although  this  term  is  much  smaller 
than  the  others. 

Although  the  preliminary  results  obtained  in  this  investiga¬ 
tion  are  promising,  and  indicate  that  it  is  possible  to  model  the 
terms  in  the  energy  equations,  and  in  particular  in  the  total  energy 
one,  accurately,  further  work  is  required  to  extend  these  results 
to  cases  in  which  the  pressure-dilatation  is  significant,  as  well  as 
to  inhomogeneous  flows.  A  posteriori  testing  of  the  models  in 
actual  calculations  is  also  necessary  for  a  complete  evaluation  of 
the  model  performance. 
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G0V6rning  eQUBtionS  E  *  ^  u,«,/2  is  tho  total  ciuiiKy  por  unit  mass,  and  the  diffusive  fluxes  :in*  by 


Here  />  is  the  density,  is  the  velocity  in  the  Xj  direction,  p  the  pressure,  f  --  r,.T  is  the  One  method  to  evaluaH*  the  p(ufor?nance  of  models  for  LES  or  RA.NS  calculations  is  the 

interna!  enernv  per  unit  mass,  T  the  temperature:  h  -  t  ♦  pj p  is  the  enthalpy  per  unit  mass:  «  priori  test,  in  wliich  the  velocity  fi(‘lds  ol)tained  from  a  <lir(‘<  t  simulation  are  Hlt<*n'd  to 


is  M'parati’ly;  This  tuixod  \va>i, also  list'd  by  ErU’hadu’r  «/’( 1992)  and  Zauj^  f.'<  «/.  ( 11)02],  aiul  is 


riiodol  of  Y'oshizawa  (1986)  introducing  the  mixed  model.  In  this  way,  tht^  eddy- viscosity  Comte  and  Lesieur  (1998)  proposed  incorporating^  this  term  into  a  modified  pressure  V. 

contribution  provides  the  dissipation  that  is  underestimated  by  purely  scale-similar  models.  This  leatls  to  the  presenct?  of  an  additional  term  in  the  equation  of  stat<?,  which  lakes  tin* 


constants  assigricd  u  priori,  sucfi  as  the  Srnaf^orinsky  (19G3)  uioflcl,  <lo  not  have  this  cliarao 
teristic).  For  internn'diate  filt(T-\vidths,  up  to  A/A  ~  IG,  tlie  l)est  prediction  is  ^iven  by  the 


5.3  SGS  turbulent  diffusion 
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since  7-(uj,  u/t)  =  Tj^.  Using  this  f«rnialisrii,  scah’-siinilar  models  can  he  derived  hy  approx¬ 
imating  the  quadratic  terms  using  the  filt<*red  velocities  Uj  to  replact^  tlse  v(^locities  uy,  for 


Tf  Si 
iii  .  H 


litiTature  to  date.  One  possibility  is  tx)  pararneteri/e  it  using  a  srah^-sirnilar model  of  the 


Multi-Block  Large-Eddy  Simulations  Proposed  niiining  head:  Multi-block  LES  of  turbulent  l)oundary  layers. 


65M06  Finite  difference  methods 


ticar-wall  n’^ion  |2j.  Tin*  highly  irn'^iilar  intrifa<<*  tin*  turbulent  anr!  noti-turbulent  iiuxieled.  LES  can  result  in  significant  savings  over  DNS,  in  terms  of  eoinputationai  costs 


How,  whifh  is  torhnoloKically  more  iiiijxirtaul,  but  <!oo.s  not  require  spt’eial  treat riietit  of 


the  euKine  nacelles  etc.).  This  issue  requires  furtfier  study  and  is  not  addressetl  here.  „,fjjis  {where  <  =  C„T  is  the  intcTiial  (uier^y  p<‘r  unit  nuuss.  C'„  is  tlu’  specific  heat  at  constant 

In  this  study  the  multi-block  approacli  is  proposed  and  tested  in  a  wall-bounded  flow  with  ^  volume  and  T  the  tcinperatiire),  p  is  the  thermodynamic  pressure,  rr^,  is  tlie  viscous  stress 

one  direction  of  homogeneity,  namely  a  temporally  developing  boundarv  layer.  This  flow  tensor  and  </j  is  the  conduction  Imat  flu.v.  In  the  pre.senf  work,  tin' fluid  is  assumed  to  be  an 

was  chosen  because  it  contains  many  of  the  features  of  importance  in  a  spatially  developing  idral  gas  with  constant  specific  lu'ats  (whose  ratio  is  7  i  ^l).  Thus,  the  pressure  is  relatecl 


turl)iiient  Mach  tunnluT  of  ().ri2.  In  their  study,  however,  the  (liveri^ence  of  ^7,  was  found  Ut  The  Kavn*-fillere<!  X a vi(^r- Stokes  equations  (10  12)  thus  heeonie 


-  (O  (“M  scheme  with  good  modified  wavenumluT  performance  and  a  small  level  of  fiissipalion  at  high 


I  S{>('()ti<)-t)r(l<‘t'ar(  uratr  fortnula.  At  tin*  t(»|»  of  the  (loiiiaitu  liir  iioriiiai  ^radirnts  of  the 
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The  first-  and  sfcond-ordcr  statistics  obtaim'd  nsin^  tin*  multi-block  approach  with  a 


for  M  =  0.3.  101,  thin  transit  ion.  (a)  0  ()1:  (1))  fuj/S*  —  0.4;  (c)  fu-,IS*  =  1.  many  times  as  necessary  in  the.  spanwise  clirc^ction.  If  nonronforminj;  im*sfjes  are?  us(*ej,  such 

Dashed  lines;  virtual  block. 

.  ^  •  that  the  spanwise  and  streaniwise  spacinj^s  of  the*  outer-layers  can  be*  increstsed  ejve*r  those?  of 


|2()|  P.  Moiu.  K  I)  S<iuirrs.  \V  H.  Caliol,  and  S.  Lrc  "A  dyiiarnic  sn])grid-scaic  inodri  for 
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C^l  ^ 


Subgrid-Scale  Models  for  Compressible  Large-Eddy  Abstract 

Simulations  pnor?  s(  ikI.v  of  sjjligri(l-sra.l<*  tjkxIcIs  for  the  n tjrl(>.sc<l  terms  in  th<* 

<’(]ti:ition  is  carried  out  iisitift  the  flow  field  ohtaituxl  from  lh<*  direc.l  sitJiiilatioti  of 
homo^enc'oiis  iso( ropie  tiirluilenee.  Seale-similar  models  iitvrdve  multiple  filtr'riiig;  op 


Tt*.  is  jluxlelfd  sej»Hr»t<d v:  I'liis  tiiixed  iiuk!<*1  was  also  used  liy  Krlebarlier  e/  nl.  (1?)!)2)  aii<l  Zaug  rl  ai  (1992).  and  is 


modelof  Yoshizawa  (1986)  iiitroduciiig  llie  mixe<j  mode).  In  tliis  way.  tlie  eddy-viseosity  (\)mte  and  Ijesleni  (1998)  proposed  iii(orpt)rat ing  Uiis  lerm  itilo  a  modified  pre.ssnre  'P. 

contril^iition  provi<Ies  tiie  <iissipation  tiial  is  underestimated  l»y  purely  scale-similar  models.  I'liis  leatls  to  the  presenre  of  an  additional  term  in  lli<*  ecpiatioii  of  state,  vvliieh  takes  the 


totisi Hilts  Hssimted  a  firion,  sui  h  as  the  SiiiHj>«>riiiskv  ( inotlel,  do  not  have  this  rliarac 


5.3  SGS  turbulent  diffusion 


'sing  this  lonnalisin,  scHle-siinilar  models  ran  lx*  derived  by  approx 


liferahire  I<j  date.  One  possihility  is  tii  pai ameterize  it  nsing;  a  srale-sitnllar  model  of  the 
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